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1 Introduction 



The first-order Riccati equation 

y' = P{x)y 2 + Q{x)y + R(x) 

is important mainly because it is a nonlinear one but directly related to the general linear 
differential equation of second-order via a Cole-Hopf transformation. It is usually considered as 
the first instance in the study of nonlinear equations pQ and is endowed with many interesting 
properties. For example, it is a Lie system admitting a nonlinear superposition principle and it 
is the only nonlinear equation of the form y' = f(x,y), where f(x,y) is a rational function of 
the variable y with coefficients analytic in x, that possesses the Painleve property (nevertheless 
the Lie-Scheffers theory or the Painleve approach will not be considered in this paper). 

The (first-order) Riccati equation is therefore a nonlinear equation that has been intensively 
studied by many authors. The important point is that it has been proved that it admits higher- 
order generalisations which are also studied by making use of several different approaches [21 [3j [4] 
(according to Davis these higher-order equations were first considered by Vessiot in 1895). All 
the higher-order Riccati equations can be linearised via a Cole-Hopf transformation to linear 
differential equations. It is known, that the higher-order Riccati equations play the role of 
Backlund transformations for integrable partial differential equations of higher-order than the 
KdV equation. The Riccati chain without potential is naturally associated to Faa di Bruno 
polynomials. The Faa di Bruno polynomials appear in several branches of mathematics and 
physics and can be introduced in several ways. 

In fact higher-order Riccati equations are related to the existence of symmetries [3 [6], Darboux 
polynomials and Jacobi multipliers (TOl ttU 12] • We also mention that the second-order 

Riccati equation has been studied in [13] from a geometric perspective and it has been proved 
to admit two alternative Lagrangian formulations, both Lagrangians being of a non-natural 
class (neither potential nor kinetic term) . An analysis of the higher-order Riccati equations and 
all these properties (Lagrangians, symmetries, Darboux polynomials and Jacobi multipliers) is 
presented in [H] . 

The Abel differential equation can be considered as the simplest nonlinear extension of the 
Riccati equation [EJ QSJ Q^. The Abel equation of the first kind [El ESI EB ED EZl E3] is given 
by 

y' = fo(x) + h(x)y + f 2 (x)y 2 + f 3 (x)y 3 . 
There is also another related equation, called Abel equation of second kind, given by 

[go(x) + gi{x)y]y' = f (x) + fi{x)y + f 2 (x)y 2 + h{x)y 3 , 

which is reducible to the previous one [231 [25] and it is not going to be considered in this paper. 
On one hand it has striking similarities with the Riccati equation but on the other side, as the 
non-linearity is of higher degree, the properties are different (and in fact more difficult to be 
studied). The objective of this paper is to study a chain of higher-order Abel equations using 
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as an approach the analysis of the differential geometric properties, the Lagrangian formalism 
and the theory of Darboux polynomials and Jacobi multipliers. 

The plan of the article is as follows: In Section 2 we review the hierarchy of higher-order Riccati 
equations and then introduce in a similar way the hierarchy of higher-order Abel equations. 
Section 3 is devoted to a particular case of second-order Abel equation. We study the existence 
of a Lagrangian formulation, obtain some constants of the motion and establish the relationship 
of this equation with the theory of Darboux polynomials and Jacobi multipliers. Section 4 
is devoted to the third and fourth-order equations and in Section 5 we consider the general 
n-dimensional case. Finally in Section 6 we make some brief comments. 



2 Riccati and Abel equations 

Let us start with the definition of higher-order Riccati equations. It is known that these equations 
can be obtained by reduction from the Matrix Riccati equation. The matrix Riccati equation 
plays an important part in the theory of linear Hamiltonian systems, the calculus of variations, 
and other related topics. 

2.1 Hierarchy of higher-order Riccati equations 

Let us denote by Dr the following differential operator, depending on a real parameter fc£R, 
that will be called 'differential operator of Riccati' 

in such a way that the action of leads to the following family of differential expressions 
B° R x = x 



Drx = ^— + kx^j x = x + kx 2 
D|x = + kx^j x = x + 3kxx + k 2 x 3 



B R x = (J^ + kx^x =x +Akxx + 6k 2 x 2 x + 3kx 2 + fc 3 x 4 
B R x = (Jjj- + kx^x = x iv) + 5kx x +Wkxx + I5k 2 xx 2 + Wk 2 x 2 x + 10A; 3 x 3 x + k A x b . 



The Riccati equation of order m of the higher-order Riccati hierarchy (o chain) , is given by 

D^x = 0, m = 0,1,2,... 

In fact, the most general form of a Riccati equation of order m is just a superposition of all the 
previous equations (linear combination the different members of the hierarchy) 

(poDl + PlB^ 1 + . . . + Pn-^R + p n )x + Pn+1 = , 
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where each pi is a function of t. 

These equations have certain properties that make them interesting from both physical and 
mathematical points of view. Next we point out some of them. 

(1) The higher-order Riccati equation of order m, member of the Riccati hierarchy, admits 
the maximal number of Lie point symmetries that can admit an equation of order m. 

(2) The higher-order Riccati equation of order m can be linearised and presented as a linear 
equation of order m + 1. 

(3) The dimensional reduction of a linear equation of order m+ 1 leads to the Riccati equation 
of order m. 



2.2 Hierarchy of higher-order Abel equations 

The most natural generalisation of the Riccati equation is 

x = f(t,x), 

where f(t, x) is a polynomial in the variable x (with coefficients depending on t). The particular 
case of f(t,x) being a cubic polynomial 

/(*, x) = Ao(t) + Ax{t)x + A 2 (t)x 2 + A 3 (t)x 3 , 

is called Abel equation. Such an equation can be considered as the simplest nonlinear extension 
of the Riccati equation. 

Let us denote by T)a the following differential operator, depending on a real parameter fceR, 
to be called 'Abel differential operator', 

B A = ± + kx 2 (t), 
at 

in such a way that the action of Da leads to a family of /c-dependent differential equations whose 
first members are given by 



T>° A x = x 



T>ax = + kx 2 ^j x = x + kx 3 

D 2 A x = + kx 2 ^j 2 x = x + Akx 2 x + k 2 x 5 

T) A x = + fcc 2 ) 3 x =x +5kx 2 x + 8kxx 2 + 9/c 2 x 4 x + k 3 x 7 

T>\x = (^- + kx 2 Yx = x iv) + 2A;(4x 3 + Uxxx + 3x 2 x) + 2/cV(22i 2 + 7xx) + Wk 3 x 6 x + k 4 x 9 



,dt 

We call this family the hierarchy of higher-order Abel equations. The Abel equation of order 
m, written in the so-called simplified form, is given by 

T>ax = 0, m = 0,1,2,... 
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Actually, the most general form of the Abel equation of order m is just a superposition of all the 
previous equations (linear combination of the different members of the hierarchy with functions 
Pi(t) as coefficients) 

{ PQ D n A + piD™" 1 + . . . + p„_iD A + Pn )x + Pn+1 = . 

An important point is that the Abel equation can not be obtained, as the Riccati equation, by 
a reduction procedure from a linear equation. 



3 Abel equation of second-order 

In this section we will analyse the particular case of the second-order Abel equation. In partic- 
ular, we describe the Lagrangian formulation of the second-order Abel equation. 



3.1 Lagrangian formalism 

The action of on the function x(t) leads to the nonlinear equation 

£ +4ta .(!) + *v-o, (i) 

that represents the Abel equation of second-order. It can be presented as a system of two 
first-order equations 

dx 

it = V 

— = —4kx 2 v — k 2 x 5 
\ at 

that determines a dynamical system that, in differential geometric terms, is represented by the 
following vector field 

T (2) =v—+F A2 — , F A2 = -4kx 2 v - k 2 x 5 . (2) 

defined on the phase space M 2 with coordinates (x,v). 
It has been proved in [13] that the second-order Riccati equation 

_ + 3 ^_) + *V = o 

can be considered as the Lagrange equation determined by the following Lagrangian 

1 



Lr 



v + k x 1 



Proposition 1 The nonlinear Abel equation of second- order HP admits a Lagrangian formula- 
tion with a non-polynomial Lagrangian. 
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Proof: There are two different ways of obtaining a Lagrangian function for the nonlinear Abel 
equation; the Helmholtz approach and the generalisation of the method used for the correspond- 
ing Riccati case. 

The Helmholtz conditions are a set of conditions that a multiplier matrix gij(x,x,t) must 
satisfy in order for a given system of second-order equations 

Xj = fj{x, x, i) , j ' = 1, 2, . . . , n, 

when written of the form 

9ijXj = gijfj(x, x, t) , i, j ' = 1, 2, . . . , n, 

to be the set of Euler-Lagrange equations for a certain Lagrangian L [261 113 US ES] (the 
summation convention on repeated indices is assumed). If a matrix solution g^ is obtained then 
it can be identified with the Hessian matrix of L, that is gij = dL/dvidvj, and a Lagrangian L 
can be obtained by direct integration of the gij functions. The two first conditions just impose 
regularity and symmetry of the matrix g^ ; the two other ones are equations introducing relations 
among the derivatives of gij and the derivatives of the functions /j. Here we only write the fourth 
set of conditions that determine the time-evolution of the g^ 

1 df a 

r(ftj) = Qik-^-kj + Qjk-Aki j A-ab = ~2~dv~ ' 

When the system is one-dimensional we have % = j = k = 1 and then the three first sets of 
conditions become trivial and the fourth one reduces to one single first-order P.D.E. 

r M+(£)«=-(£)+'(lM§£)'-°- « 3 > 

that in the case of the Abel equation becomes 

_ {4 kx 2 v + k 2 x 5 ) (^) - Akx 2 g = . (4) 

So, the problem reduces to find the function g as a solution of this equation. Once a solution 
g is known a Lagrangian L is obtained by integrating two times the function g. The funtion L 
obtained from g is unique up to addition of a gauge term 

Next we consider the second method that is specific for this particular nonlinear problem. 
The starting point is the idea that, since the Abel equation is very close related with the 
Riccati equation, it seems natural to assume that the Abel Lagrangian must be a non-polynomial 
function similar to that of the second-order Riccati equation. 

Let us begin by considering the following one degree of freedom Lagrangian 

L = (v + kU(x,t)) m ' (5) 
From such a Lagrangian we arrive to the following second-order nonlinear equation 

^LLH!l \ kU x x + ( — - — ) k 2 UU x + kll t = . (6) 
1 + mJ \l + mJ x 1 y ' 
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Hence, in the particular case of U and m being given by 



U(x,t) = x', m = 2, 

then the Lagrangian ([5]) leads to ([T]). Thus the second-order Abel equation ([1]) turns out to be 
the Euler-Lagrange equation of the Lagrangian function 

LA= (v + lxY' (7) 

Finally, as a byproduct of this approach, we have also obtained the Lagrangians for the whole 
family of nonlinear equations ((6j) depending of a function U. □ 

As a corollary of this proposition we can state that when the function U is time-independent 
the nonlinear equation © has a first-integral that can be interpreted as a preserved energy. That 
is, if we restrict the study to nonlinear equations arising from a time-independent Lagrangian 
of the form 

L 



(v + kU{x)) m 

then we can define an associated Lagrangian energy El by the usual procedure 

E L = A(L) -L, A = v , 

and we arrive to 

_ - ((l + m)v+ kU(x)) d_ 
L (v + kU{x)) m + l ' dt L 

In the particular case of the Abel Lagrangian La we have 

_ (3v + kx 3 ) d _ 

Ela ~- ( V +kx*r> jt ELA ~°- (8) 

Note that La is non-natural and, as there is neither kinetic term T nor potential function V, 
the energy cannot be of the standard form El = T + V. But, in spite of its rather peculiar 
form, El a is a conserved function for the Abel equation. 

An important property of the Lagrangian formalism is that for one degree of freedom systems 
if an equation admits a Lagrangian formulation then the Lagrangian is not unique [37^138]. This 
property can be proved in two different ways. First, the Helmholtz equation §S§ is a linear 
equation in partial derivatives and thus it admits many different particular solutions. Moreover 
it is clear from the form of the equation © that if g\ is a particular solution then g<i = fg\ with 
r(/) = is also a solution. A second method is related with the properties of the symplectic 
formalism. In a two-dimensional manifold all the symplectic forms must be proportional. Hence 
if ul is known then any other symplectic form u)2 must be proportional to ul, that is u)2 = f^L- 
Then 

i(T L ) lo 2 = f i(T L ) uj l = f dE L . 

The right-hand side is an exact one-form if, and only if, dfAdEL = 0, which shows that / must 
be a function of El- In this case it can be proved that the new symplectic form u>2 is derivable 
from an alternative Lagrangian L^L for Tl- 
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In the particular case of the Abel system T^ 2 \ several alternatives Lagrangians can be obtained 
that, in most of cases, are of non-algebraic character (with logarithmic terms). Nevertheless, in 
the particular case of / given by / = (— l/£x A ) 4//3 , we have obtained the following algebraic 
function 

L A = (3v + kx 3 ) 2 / 3 (9) 

as a new alternative Lagrangian for the Abel equation ([1]). This new Lagrangian is equivalent 
to La in the sense that both determine the same dynamics. It determines a new energy E^ A 
that is a constant of the motion for the Abel equation; nevertheless it must not be considered 
as a new fundamental constant since it is a function of the original energy E^ A . 

3.2 Constants of the motion and geometric formalism 

A function T that satisfies the following property 

dt r dt m r dt m+1 

is called a generator of integrals of motion of degree m. Notice that this means that the function 
T is a non-constant function generating a constant of motion by successive time derivations. 

(2) 

Let us denote by T\ the following function 



T 



(2) 



v + k x 3 

Then we have that under the evolution given by Abel's equation ([T|) 

! T (2) =T (2) A r f = . 

dt 1 2 dt 2 



Thus, the function J t \ defined by 



Jti — 1 1 —t. 



is a time-dependent constant of the motion for the Abel equation. 

This means that we have obtained two constants of the motion (of quite different nature) for 
the Abel equation of second-order: the energy El a and the time-dependent function Jt\- 

In differential geometric terms a time-independent Lagrangian function L determines an exact 
two- form u>l defined as 

Ol = yQ~) dx ' ul = - d6 L , 

and L is said to be regular when the 2-form ujl is symplectic. In the particular case of L given 
by © oj La is given by 
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and the dynamical vector field I^ 2 ) is the solution of the equation 

i(r®)u LA = dE LA . 

Next we consider two interesting class of symmetries: 'master symmetries' and 'non-Cartan 
symmetries'. The idea is that, in differential geometric terms, constants of motion that depend 
of the time but in a polynomial way are related with the existence of master symmetries [3Q|, 
[3lj [32j [35] and in some very particular cases with non-Cartan symmetries. 

Given a dynamics represented by a certain vector field T, then a vector field Z satisfying 

[Z,T]=Z^0, [Z,T]=0, 

is called a 'master symmetry' of degree m = 1 for T. When Z is such that 

[Z,T]=Z^0, [Z,T]^0 and [[Z,T},T}=0, 

then Z is called a 'master symmetry' of degree m = 2. The generalisation to higher values of m 
is straightforward: 

(ad(r)) m+1 (Z) =0, but (ad(r)) m (Z) ^0. 

It is well-known that symmetries are important because they give rise to constants of the 
motion and reduction procedures. Master symmetries, which are a rather peculiar class of 
symmetries, determine time-dependent constants of motion (the system is time-independent 
but the constant is however time-dependent). This can be seen as follows: if Z is a master 
symmetry of degree one, the time-dependent vector field Yz determined by Z as follows |33j 

Y z = Z + t[Z,T] + ^)t 2 [[Z,T],T] 

is a time-dependent symmetry of T t = d/dt + T, which is the suspension of the vector field 
r [33]. This symmetry determines a time-dependent constant of motion Jt = T — tT(T) that 
depends linearly of t (for m = 2 the corresponding constant Jt will be quadratic in t and for 
m = 3 will be cubic). 

(2) 

Let Z\ be the Hamiltonian vector field of the function T± , that is, the unique solution of the 
equation 

i{Z 1 )u JL = dT[ 2) , 

which is given by 

Z 1 = -(l)Pl 1 (x-^ + (v-2kx s )-^), P A1 = v + kx\ 
Then Z\ is a symplectic symmetry (that is, Lz x = 0) because it is the Hamiltonian vector 

(2) 

field of T-y , and moreover it is a dynamical symmetry because 

, T^]lo l = i(Zi)(£ r(2) u; L ) - £ r(2) (i(ZiK) = -£ r(2) (dTf >) = 0, 
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and therefore, as ujl is non-degenerate, \Z\ , r^ 2 )] = 0. 

Note however that Z\ is not a symmetry of the energy since Z\(El a ) ^ 0. 

Thus, Z\ is a dynamical but non-Cartan symmetry of the Lagrangian system [35} 136] , These 
symmetries are rather peculiar and only appear in some very particular cases. In particular it 
was proved in [36] that if the Hamiltonian vector field Xp with the function F as Hamiltonian 
in a symplectic manifold (M, lo) is a dynamical but non-Cartan symmetry, then Xp(H) must 

(2) 

be a numerical constant Xp(H) = a/ 0. In this case we are considering, F = T} ' and we have 
Z 1 (E hA ) = a = -\. 

We close this section by recalling that the Riccati equation was endowed with similar properties 
but the function T x was the Lagrangian Lr itself [H] . 

3.3 Darboux polynomial and Jacobi multiplier approach 

The existence of constants of the motion and the Lagrangian inverse problem for polynomial 
vector fields are two questions related with two important ideas: Jacobi multipliers and Darboux 
polynomials. 

Let U be an open subset of W 1 . We say that a polynomial function T> : U — > R is a Darboux 
polynomial for a polynomial vector field X if there is a polynomial function / defined in U such 
that XV = fV [H [9j [H]. The function / is said to be the cofactor corresponding to such a 
Darboux polynomial and the pair (f,T>) is called a Darboux pair. 

When / = 0, then the Darboux polynomial is a first integral. We say that T> is a proper 
Darboux polynomial if / 7^ 0. If T>\ and T>2 are Darboux polynomials with the same cofactor, 
the quotient T>x/T>2 is a first integral. 

On the other side given a vector field X in an oriented manifold (M,Q), a function R such 
that Ri(X)Q is closed is said to be a Jacobi multiplier (JM) for X. Recall that the divergence 
of the vector field X (with respect to the volume form f2) is defined by the relation 

C x Fl = (divX)O. 

This means that R is a multiplier if and only if R X is a divergenceless vector field and then 

c RX n = (div rx) n = [x(R) + R dwx] n = o, 

and therefore we see that R is a last multiplier for X if and only if 

X(R) +RdivX = 0. (10) 

Note that if R is a never vanishing Jacobi multiplier, then fR is a Jacobi multiplier too if and 
only if / is a constant of motion. We also note that the equation ([4j) can now be considered as 
a particular case of the equation (I10p . 
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The remarkable point is that if T>i, . . . ,T>k, are Darboux polynomials with corresponding 
cof actors fi, i = 1, . . . , k, one can look for multiplier factors of the form 



R = f[v? (11) 



i=i 

and then 

i=i 1 i=i 

and therefore, if the coefficients Vi can be chosen such that 

k 

Y,"ifi = -dwX (12) 

i=l 

holds, then we arrive to 

i=l 

and consequently i? is a Jacobi last multiplier for X. 

Finally, if -R is a Jacobi multiplier for a vector field which corresponds to a second-order 
differential equation, there is an essentially unique Lagrangian L (up to addition of a gauge 
term) such that R = d 2 L/dv 2 [El QU [Lj . 

From these general concepts we can return to the Abel equation. In this case the polynomial 
T>\ defined by 

T>i(x, v) = v + kx 3 
is a Darboux polynomial for r*- 2 -* with cofactor —kx 2 since 

d d \ 

v — + f A2tt) ( v + kxS ) = ~kx 2 (v + kx 3 ) . 
ox ov J 

The divergence of the vector field is —4kx 2 , and then, according to (|12p . we see that there 
is a multiplier of the form 

R = V\ X 



with v\ = —4. Consequently, the Abel equation admits a Lagrangian description by means of a 
function L\ such that 

d 2 Li 34 

-g^r = ( v + kx ) ' 

from where we obtain the Lagrangian Li = given by ([7]). 
But the polynomial T>2 defined by 

T>2(x, v) = 3v + kx 3 
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is a Darboux polynomial for I^ 2 ) with cofactor —3kx 2 , because 
d d \ 

V dx + FA2 dv ) ^ + = 3kx2v ~ 3 ( Akx2v + ^ x > = ~3 kx2 ( 3v + kx3 ) i 

and then, using the equation (|12p we can find another Jacobi multiplier of the form T> v 2 2 with 
V2 = —4/3. The Abel equation admits a Lagrangian description by means of a function L2 such 
that 



d 2 u 



2 =(3v + kx 3 ^ 4 / 3 



dv 2 

from where we obtain the Lagrangian L2 = La given by ([9]). 

Remark that, as indicated above, if P and Q are two Darboux polynomials with the same 
cofactor then Pj Q is a constant of the motion. This is just what happens with the energy El a 
obtained in (jHJ) which is given by V 2 /V\ (up to the sign J . 



4 Abel equations of third and fourth-order 

In the following all the functions that appear as constants of the motion will depend on the time 
t. So we better consider to describe them as first integrals that depend polynomially on time. 



4.1 Abel equation of third-order 



The action of the operator T)a three times on the function x{t) leads to the following nonlinear 
equation 



d 3 x r > /d"x\ 

1_ Zkr 2 



dt 3 



5kx 



\dt 2 J 



+ 8kx(^-) +9k 2 x 
\ dt J 



2 .4^^ 



+ k 3 x 7 



0. 



(13) 



that represents the third-order element of the Abel equation chain. It can be presented as a 
system of three first-order equations 



dx 

v 

t 
a 

~dl 



(14) 



-5kx a 



i-XV 



9k 2 x 4 v — k 6 x 



3^7 



that represents a dynamical system that, in differential geometric terms, is represented by the 
following vector field in the phase space M 3 , with coordinates (x,v,a) 



d d „ d 

r = v ir + a ir + F ^ir- 

ox ov oa 



F A3 = -5kx 2 a - 8kxv 2 - 9k 2 x 4 v - k 3 x 7 . 



(15) 



In what follows we make use of the following polynomials 



Pao = x . 



2^.5 



Pax = v — Fai = v + k x , P A 2 = a — Fa2 = a + ^kx v + k x 
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defined on the phase space and obtained by making use of the substitution x \— > v and Jho. 
Then we have 

r (3) (-FUo) + kx 2 P A0 = v + kx 3 , 

T ( - 3) (P A1 ) + kx 2 P M = a + 4kx 2 v + k 2 x 5 , 

T {3) (Pa2) + kx 2 P A 2 = F A 3 + 5kx 2 a + 8kxv 2 + 9k 2 x*v + k 3 x 7 . 

that can be rewritten as follows 

Y^{P A0 ) + kx 2 P A0 = Pai, 

Y^{P A1 ) + kx 2 P A1 = P A2 , 

T i3) {P A2 ) + kx 2 P A2 = 0. 

Note that according to these properties P A2 is a Darboux polynomial with / = —kx 2 as cofactor. 
The divergence of the vector field is — 5kx 2 , and using relation (|12p we see that R = (P A2 )^ 2 
with H2 = — 5 is a Jacobi multiplier. 

(3) 

Next let Tj be the following function 

T (3) _ x 
1 ~ PA2' 

and then we have 

r<s>(zf >) = rf = v + kx " , r( 3 )(rf) = t 3 (3) = i , r( 3 )(T 3 (3) ) = rf = o . 

PA2 

This means that and T 2 ^ are generators of of constants of motion for the third-order 
element of the Abel equation chain represented by the dynamical vector field T^ 3 \ Thus we can 
state the following proposition. 

Proposition 2 The two functions Jti and Jt 2 defined as 

j tl = r 2 (3) -t, j t2 = if - 1 r 2 (3) + {\)t\ 

are first integrals, that depend polynomially on time, for the Abel equation of third- order. 

Note that Jti is linear in the time t and Jt 2 is quadratic. So these expressions are similar to 
the constants of the motion determined by master symmetries; nevertheless in this third-order 
case we have not made use of any symplectic structure and we have obtained these functions 
without relating them with symmetries of a symplectic structure. This is an interesting situation 
deserving an additional analysis in the next sections. 

Note also that both, Jti and Jt 2 , can be written as quotients of polynomials; so if we consider 
the system as a time-dependent system then the dynamics is geometrically represented by the 
vector field = + d/dt and the following polynomials 

V 2 = P A i - tP A2 , V 3 = P M - tP A1 + {^)t 2 P A2 , 
are two Darboux polynomials with the same cofactor as P A2 

rf (Vi) = (r< 3 ) + ^) (pi) = - kx 2 v u i = 2, 3. 
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4.2 Abel equation of fourth-order 

The action of four times on the function x(t) leads to the following nonlinear equation 

x iv) + 2A;(4x 3 + Uxxx + 3x 2 x) + 2k 2 x 3 (22x 2 + 7xx) + 16k 3 x 6 x + k 4 x g = , (16) 

that represents the fourth-order element of the Abel equation chain. This equation determines 
a dynamical system that, in geometric terms, can be represented by the following vector field 
on M 4 as phase space, with coordinates (x,v,a,w): 

r (4) =^ + a# + ^ + i^, (17) 
ox ov oa aw 

where 

F A4 = -2k(Av 3 + Uxva + 3x 2 w) - 2k 2 x 3 {22v 2 + 7xa) - Wk 3 x 6 v - k 4 x 9 . 

Now we introduce the polynomial P A3 

P A3 = w - F A3 = w + 5kx 2 a + 8kxv 2 + 9k 2 x 4 v + k 3 x 7 , 

obtained from the expression of x with the substitution x \— > v, x > a y i h «j. Then we 
have 

T^(P A0 ) + kx 2 P A0 = v + kx 3 

T^{P Al ) + kx 2 P Al = a + 4kx 2 v + k 2 x 5 

T {4) {P A2 ) + kx 2 P A2 = w + 5kx 2 a + 8kxv 2 + 9k 2 x 4 v + k 3 x 7 

T (4 \P A3 ) + kx 2 P A3 = F A4 + k{8v 3 + 2Qxva + Qx 2 w) + k 2 (Ux 3 v 2 + 14x 4 a) 

+ 16k 3 x 6 v + k 4 x 9 



that can be rewritten as follows 



T^(P M ) + kx 2 P M = P A1 

T^(P Al ) + kx 2 P Al = P A2 

r^(p A2 ) + kx 2 p A2 = p A3 

T^(P A3 ) + kx 2 P A3 = 



Let now 

T (4) 

be the following function 



(4) 



1 " Paz 
and then we have 

r (4 )(T (4 )) = T (4) rWpW) = T w r W(r 3 (4) ) = t 4 (4) , rW(r 4 (4) ) = o , 

with , Tg 4 "* , and given by 

(4) _ v + kx 3 (4) _ a + 4kx 2 v + k 2 x 5 (4) _ w + . ..+ & 3 ^ 7 _ 
-PA3 3 Pa3 4 -P.A3 
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Proposition 3 The three functions Jti, Jt2, o,nd Jt 3 defined as 

j t i = n-t 

Jt2 = T 2 -tT 3 + (^)t 2 

J t3 = Tl -tT 2 + (^)t 2 T 3 -(^)t 3 

are first integrals, that depend polynomially on time, for the fourth-order element of the Abel 
equation chain. 

The situation is similar to the n = 3 case and the functions Jt r , r = 1, 2, 3, are polynomials of 
order r in the variable t. 



5 Equation of Abel of order n 



We have seen that the second-order element of the Abel equation chain is endowed with some 
specific properties (e.g., it admits a Lagrangian description) but the third and fourth-order 
elements of the chain also enjoy very similar properties. Now in this section we study the 
equation of order n and prove that these properties characterise to all the equations of the 
family in an independent of the order way. 

The equation of Abel of order n can be obtained as the equation arising from the action of 
the operator on the equation of order n — 1 

D A (D^ _1 x) =D^x = 0. 

This equation determines a dynamical system that, in geometric terms, can be represented by the 
following vector field defined on the phase space M n , with coordinates (x = x±, x 2 , x 3 , . . . , x n ): 

^(n\ d d d „ d , . 

r (n) =* 2 — + x 3 — + x A — + ... + F An — , 18 

OX OX2 ux 3 ox n 

where F^n is obtained a from the expression for D^x with the substitution x 1— > x\, x — ► x 2 , 

In the previous sections we have made use of the polynomials Pao, Pai, Pa2, and P43 defined in 
the phase space and whose explicit expressions, when written in the notation of the coordinates 
xi,x 2 ,x 3 , ... ,x n , were given by 

Pao = x x , 

Pai = x 2 - F A i = x 2 + kx\, 

Pa2 = x 3 - Fa2 = x 3 + Akxfx 2 + k 2 x\ , 

P A3 = X4 — Fa 3 = X4 + 5kx 2 x 3 + 8kx\x 2 + 9k 2 xfx 2 + k 3 x\ . 

In the general case we have P A n-i = x n — F An -i that leads to an expression of the form 

P An -\ = x n - F An -i = x n + (n + l)kx 2 x n - 1 + ... + k n ~ 1 x 2n ^ 1 . 
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Proposition 4 The action of the dynamical vector field on the n polynomials P Ar , r = 
0, 1, 2, . . . , n — 1, is given by 

(i) T^\P Ar ) + kx 2 P Ar = P Ar+1 , r = 0,1,2,..., n- 2, 

(a) rW(P An „ 1 ) + kx 2 p An - 1 = o. 

The property (i) can be proved by induction. The property (ii) follows by direct calculus. □ 
Proposition 5 The time evolution of the functions P Ar /P An -\ is given by 

rK^)=i^±l, r = 0,l,2,...,n-2. 

By direct calculus we have 

M ( PAr \ = T {n) {P Ar )P A n-l - P Ar T {n) (PAn-l) 

Then, making use of the properties (i) and (ii) of proposition (jj]), we arrive at 

r (n) ( P Ar \ _ [PAr+l ~ kx 2 P Ar ]P An -l ~ P A r[~kx 2 P An -i\ _ P Ar+ l 
\PAn-J~ {PAn-l? ~ PAn-l' 

Notice that the first and the last derivatives in this series, corresponding to r = and r = n—2, 
become 

d / x \ Pai d / P An -2 \ 

dt V P An - 1 / Pau-X ' dt V P An -i ) 

(n) 

Now let I} be the following function defined by 

r (n) 



1. 



1 PAn-l ' 

and then, making use of the two preceding Propositions we can obtain the values of the sequence 

(n) 

of time derivatives of the functions Tt , which are given by 



p(n)^y( n )^ _ >-p 



( n ) _ x 2 + kx 3 P A 



PAn-l PAn-l 



{n) („) ( n ) _ x 3 + 4kx 2 x 2 + k 2 x b _ P A2 

1 l J 2 ) — J 3 — p, — p, 

rAn-1 *An-l 

(„), („), („) _ x 4 + 5kx 2 x 3 + 8fcxx^ + 9fc 2 x 4 x 2 + k 3 x 7 P A3 

1 U3 ) — 1 4 — p~ — p 



rW(TW) = o. 



From here we can state the existence of a family of n — 1 first integrals depending polynomially 
on time. 
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Proposition 6 The (n — 1) functions Jt r , r = 1, 2, 3, . . . , n — 1, defined as the following poly- 
nomials of order r in the variable t 



■V 

:^ 2 e\-(^ 3 



Jtl 


_ rp{n) 

— 1 n-l 


- t 


Jt2 


rp(n) 

— 1 n-2 


.rp(n) 

ll n-l 


Jt3 


_ rp(Tl) 

— 1 n-3 


+ rn{n) 
ll n-2 


'tn-1 


_ rp( n ) 


-tr 2 (n) - 



:^t 2 Tt ] - + (-in^)t n 

are n — 1 functionally independent first integrals, that depend polynomially on time, for the Abel 
equation of order n. 

An alternative form of proving the existence of all these constant of the motion is as follows. 
The n polynomials T> a , a = 1, 2, . . . , n, defined in the extended phase space R™ x R as 

= Pau-i , 

T^2 = PAn-2 — tPAn-1 , 

£>3 = PAn-3 - tPAn-2 + (~)t 2 PAn-1 , 

T^i = PAn-A — tPAn-3 + (t:)* 2 PAn-2 ~ (77) t 3 PAn-1 , 

2 b 



V n = P A0 -tP Al + ... + {-l) n (^)t n P An „ l 



rl n \v a ) = frW + p, (v a ) = -kx 2 v a , a = 1, 2, 



are n Darboux polynomials with the same cofactor 

„, kx 2 V a , a = 1,2,... 

at J 

Hence the functions 

P>a 

Jtab = ^r, a,b = 1,2, . . . ,n, 

are constants of the motion. In fact, we can arrange all these functions as the entries of an n- 
dimensional matrix [ Jtab] that becomes a matrix formed by constants of the motion (the diagonal 
elements are just ones) with the fundamental set of functions Jtk placed in the first row. 

Finally, the divergence of the vector field is given by divr( n ) = — (n + 2)kx 2 . Thus, using 
relation (|12p . we obtain the following Jacobi multipliers for the Abel equation of order n (or for 
the dynamical vector field r( n )) 

R a = {VaY n , Mn = -(n + 2), a = l,2,...,n. 

We note that all these n Jacobi multipliers are different, that is Rb ^ R a , b ^ a, but they are 
proportional by a function that is an integral of the motion. 
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6 Final comments 



We have studied a chain of higher-order nonlinear Abel equations using, as starting point, the 
idea that they have many similarities with the higher-order nonlinear Riccati equations. We 
have made use of the Lagrangian formalism (inverse problem, non-polynomial Lagrangians, 
nonstandard symmetries) in the case of the second-order equation and of other mathematical 
tools (Darboux polynomials and Jacobi multipliers) in the case of higher-order nonlinearities. 
All these questions seems to be really interesting and we think they deserve a deeper study. 

Finally, we mention that all these equations have (for any order of the equation) a family of 
first integrals Jtk that depend of the time as a polynomial in t. In the symplectic case functions 
of such a class are associated to master symmetries of the (Lagrangian or Hamiltonian) system, 
but in the general Abel case we have proved the existence of such constants without refering to 
any symplectic structure. This is in fact a very interesting fact that is to be studied. 
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